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Колесные пары, непосредственно воспринимающие силовое 
воздействие от взаимодействия колес с рельсами, являются 
интенсивными источниками шума при движении транспорт- 
ных средств, например, мостовых кранов и железнодорожного 
подвижного состава. Несмотря на различия конструктивного 
исполнения и функционального назначения, расчет виброаку- 
стических характеристик может быть выполнен с единых по- 
зиций. Следует отметить, что точность расчета уровней звуко- 
вого давления и вибраций в значительной степени определя- 
ется точностью задания частотно-зависимых коэффициентов 
потерь колебательной энергии. Теоретические исследования 
этой физико-механической величины представляет значитель- 
ные трудности. Поэтому для инженерных расчетов уровней 
звукового давления и вибраций целесообразно использовать 
регрессионные зависимости, полученные на основе математи- 
ческой обработки экспериментальных величин коэффициен- 
тов потерь узла колесной пары и аппроксимация зависимостей 
степенными полиномами различных степеней. 


Ключевые слова: коэффициенты потерь, энергия колебаний, 
колесные пары. 
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Введение. Исследованию шума и вибрации подвижного состава посвящены работы [1—11]. Расчет вибраций 


узлов колесных пар технологических машин различного функционального назначения (электровозы, тепловозы, мо- 


стовые краны и т. д.) основан на использовании методов энергетического баланса [12]. Для этого необходимо знать 


значения коэффициентов потерь колебательной энергии, являющимися частотно-зависимыми величинами. В лабора- 


торных условиях проводилась серия экспериментов по определению коэффициентов потерь колебательной энергии 


узлов колесных пар. Схемы установки данных узлов представлены на рис. 1. Колесная пара устанавливалась в центрах 


на массивной плите. 
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Рис. 1. Схема установки для определения коэффициентов потерь колебательной энергии колесной пары: 
а) колеса; 6) оси колесной пары: 
1 — генератор сигналов возбуждения; 2 — анализатор; 3 — вычислительный комплекс; 4 и 5 — усилители сигналов; 
6 — исследуемый элемент; 7 — вибродатчик; 8 — виброизолирующие стойки; 9 — возбудитель колебаний 


Результаты измерений, представленные на рис. 2, показали, что при общей закономерности характера частот- 
но-зависимых коэффициентов потерь колебательной энергии их величины для узлов колесных пар существенно отли- 
чаются от известных данных для пластинчато-стержневых конструкций [13]. Поэтому полученные результаты пред- 
ставляют собой интерес для расчетов спектров вибрации и шума основных звукоизлучающих элементов общей аку- 


стической системы узлов колесных пар. 
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Рис. 2. Коэффициенты потерь колебательной энергии узлов колесных пар и режущего инструмента 


Для обоснования рациональных параметров конструкции системы шумозащиты рассчитывался количествен- 
ный вклад основных источников шума (обрабатываемых колес, оси колесной пары) в формирование спектров шума в 
рабочей зоне оператора. Как показали теоретические исследования, точность расчета уровней звукового давления со- 
ставляющих спектра в значительной мере зависит от величин частотно-зависимых коэффициентов потерь колебатель- 
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ной энергии. Для этих целей по экспериментальным данным, приведенным на рис. 3, получены регрессионные зави- 
симости коэффициентов потерь колебательной энергии колес, оси колесных пар. 

Регрессионный анализ позволяет аппроксимировать кривые, построенные по экспериментальным данным, 
или таблично заданные функции методом наименьших квадратов. Метод наименьших квадратов позволяет получить 
наиболее точные значения параметров эмпирической зависимости и, в соответствии с этим, максимально приблизить 
теоретические результаты расчетов к фактическим. 
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Рис. 3. Частотная зависимость коэффициентов потерь колебательной энергии 


Как правило, этот метод используется для нахождения линейной зависимости. Однако, метод наименьших 
квадратов и процедура подбора прямой регрессии полностью переносятся и на случай нахождения криволинейной 
зависимости при условии, когда уравнение кривой может быть после некоторых преобразований сведено к линейному 
тренду. В практике криволинейного выравнивания широко распространены два вида преобразований: натуральный 
логарифм (№) и обратное преобразование (1/х). При этом возможно преобразование как зависимой переменной у, так 
и независимой х или одновременно и той, и другой (Таблица 1). 

Таблица 1 
Кривые, сводящиеся к уравнению прямой преобразованием исходных данных 


Название кривой Уравнение Преобразование 
Экспоненциальная у = ае“" У =Шу 


Степенная у=ах” 
а 


Гиперболическая [ типа У=а+-— 
х 


Гиперболическая П типа 


Гиперболическая Ш типа 


Логарифмическая 


5-образная 


Обратнологарифмическая 





Вестник Донского государственного технического университета 2016, 1(84), 127-135 





Кроме этого, метод наименьших квадратов позволяет определить эмпирическую зависимость и в виде поли- 
нома (многочлена) любой степени т. При выборе наилучшей кривой оценивается такой параметр, как среднеквадра- 
тичное отклонение (СКО), который показывает степень приближения эмпирической кривой к исходной и поэтому 
должен быть минимальным: 

Га | 
с„= У (в, (м)-,), (1 





где п — количество исходных данных; т — степень аппроксимирующего полинома; Р.(х;) — значение эмпирической 


функции; у; — значение исходной функции. 
2 

Так как многочлен Р, (х) =а+алх+ах“ +..+а»„х” определяется своими коэффициентами, то фактически 

нужно подобрать набор коэффициентов а, а1.....а» ‚ минимизирующий функцию 

2 
ГД) И т . 
ы 2 _ ] 

Ф(ау,а....а„)= (в, у) = >. | А Е & (2) 


ОФ 
Используя необходимое условие экстремума — =0,К = 0,1,...т получаем нормальную систему метода 


бах 
наименьших квадратов: 
т Г | Ш 
ХУ" И, = Ху, К=01,.. т. (3) 
7=0\ 1=0 1=0 


Полученная система есть система алгебраических уравнений относительно неизвестных а1, а›,...‚ав, решение 
которой даст коэффициенты эмпирических зависимостей. 

При аппроксимации таблично заданной функции полиномом его степень заранее неизвестна, поэтому возни- 
кает проблема выбора оптимальной степени аппроксимирующего многочлена. Вычислим для каждого полинома ве- 


личину |“. . За оптимальное значение степени полинома предлагается принять то значение т, начиная с которого 


величина а» стабилизируется или начинает возрастать. 


Для определения аналитических зависимостей коэффициентов потерь от частоты необходимо аппроксимиро- 
вать каждую кривую отдельно. Проведем расчет коэффициентов уравнений регрессии для всех вышеперечисленных 
функций последовательно для оси колесной пары и колеса. 

Ось колесной пары. Исходные данные коэффициента потерь колебательной энергии оси колесной пары 
представлены на рис. 3. Регрессионный анализ проводим аналогично фрезе, результаты сводим в табл. 2, а также 
представляем на рис. 4, 5. 

Таблица 2 
Зависимости коэффициента потерь колебательной энергии для оси колесной пары 


Название кривой Уравнение 


0 
0 
Гиперболическая [ типа т 3,086Е-—0 
| В б П = И. ИИ 6,007Е—0 
ипероолическая ИП типа м 234.79+ 0.0061 Е : — 
0 


Гиперболическая Ш типа | П= т. 3.814Е— 


Логарифмическая |= 0,0061 — 0,00031-ш(Х) | 3,461Е-0 
, 0 


5-образная _ г 3,389Е— 
Об ЕЕ Е 3,066Е-0 
ратнологарифмическая | П 143.07 +16.51. г 1) Е 


4 
4 
4 
4 
4 
4 
4 
4 





Бондаренко В.А. и др. Экспериментальные исследования коэффициентов потерь энергии колебаний 


И 
0,0055 


—— ось 


\ 


___ Экспонента 


0,005 —_ Степенная кривая 


ПА 
"|о 


7 


00045 —*— Гиперболическая кривая | типа 


< 


Г 


У | 


—— Гиперболическая кривая || типа 


0,004 
——_ 5-образная кривая 


0,0035 


—+— Обратнологарифмическая 
кривая 


—=— Логарифмическая кривая 


АА 
ЕЕ^ 


и 





0,003 
31,5 63 125 250 500 1000 2000 4000 8000 ——-  Гиперболическая кривая Ш типа 


Рис. 4. Графики коэффициента потерь колебательной энергии оси колесной пары и аппроксимирующих кривых 
Таблица 3 


Аппроксимация коэффициента потерь колебательной энергии для оси колесной пары 
полиномами различных степеней 


ско (а, 
п= 0,01-0,004( г)+0,0005(е 2} 
п=0.005+0,003( /)-0,0018(® г}? +0,00029(в /}} 

4 п=-0,18+0,04е /)-—0,024(е Г} +0,006(е г} —0,0005е у}! 

‚66Е-05 

05 
5 


п=0,047 0,095 (в /)+0,084(в /} — 0,0360 +0.008 (2 \- | т 
— 0,0006 Г} 
П=-0.035+0,108(в /)-0/121(= л} +0,072(5 У). — 0.0247 | зе 
+0,004(= г}? —0,0003(® у} | 
т 


п=-0,301+0,875(в /)--1,0510 1) +0,686(= 1) — 0,2636 77+ | ео 
+0,059(= /}’-0,007(в д} +0,0004( г); | 
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Рис. 5. Результаты аппроксимации коэффициента потерь оси колесной пары полиномами 
Для коэффициента потерь колебательной энергии оси колесной пары по значениям СКО наиболее прибли- 
женным является полином шестой степени: 
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1=-0,035+0,108( /)-0,121(= ху +0,072(® г)’ -0,024( г}! +0,004(® г}? —0,0003( г}. (4) 
Колесо. Исходные данные коэффициента потерь колебательной энергии колеса представлены на рис. 3. Регрессион- 
ный анализ проводим аналогичным образом, результаты сводим в табл. 4 и 5, а также представляем на рис. 6, 7. 


Таблица 4 


Зависимости коэффициента потерь колебательной энергии для колеса 


н_30вл0зе 35" [59428-04 
=0.0063у и | 27338.04 


в 
Гиперболическая [1 типа И = 0,0034 + —— 2.623Е-04 
Гиперболическая П типа — 554 94400098 54.94+ 0.0091 5,972Е-04 


Гиперболическая Ш типа —577 —3775.23+297.49 7 5 г 297 49 Я 3,671Е-04 


15,11 
о 
Об = 2,5 1Е-04 
ратнологарифмическая | 771-2142. (Е й 


—Ф— колесо 
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Рис. 6. Графики коэффициента потерь колебательной энергии колеса и аппроксимирующих кривых 


Таблица 5 


Аппроксимация коэффициента потерь колебательной энергии для колеса полиномами различных степеней 


|= 0,009 —0,004(® /)+0,0005( у} 1,20Е-04 


|= 0,007 +0,0004(® /)—0,0007( у’ +0,0002(в г} 1,11Е-04 
п=-0,01+0,028( /)-0,017( Х} +0,004( у} —0,0004(= 7} | 1,14Е-05 
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п=-0,011+0,029( Г)-0,019(= Х}’ +0,005( г}’-0,0005( Г} + 
+0,00001(е Г} 
п=-0,042 +0,106(е /)- 0,096 Г} +0,046(в у} —0,012(в ду + 


1,31Е-05 


7,88Е-06 
+0,002(® } —0,0001(е ХУ 


п=-0,052 + 0,136( г)-0,133(® г)’ + 0,07(е Х}—0,022(= р} + 
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Рис. 7. Результаты аппроксимации коэффициента потерь колеса полиномами 


Для коэффициента потерь колебательной энергии колеса по значениям СКО наиболее приближенным являет- 
ся полином шестой степени: 


п=-0,042+0,106(® г)-0,096(® Г’ +0,046(в г)’ -0,012(® г}! +0,002(= г} —0,0001(® (5) 


Расчетные спектры шума представлены на рис. 8. 
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Рис. 8. Расчетные спектры шума 


Выводы. Излучение звука колес в широком частотном диапазоне 500-8000 Гц превышает санитарные нормы 
на 3—10 дБ. В седьмой, восьмой и девятой октавах превышение санитарных норм уровня звукового давления состав- 
ляет 8—10 дБ, что на 2—3 дБ больше, чем в случае оси колесной пары. Разница расчетных и экспериментальных данных 
по величинам превышения уровней звукового давления над предельно допустимыми значениями составляет 2—2,5 ДБ, 
что является вполне приемлемым для технических расчетов уровней шума и вибрации. 
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